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ABSTRACT 

The effect of liquid surface motion on the vapor condensation 
in a tank mixed by an axial turbulent jet is numerically 
investigated. The average value (over the interface area) of 
the root-mean-squared (rms) turbulent velocity at the inter- 
face is shown to" be linearly increasing with decreasing 
liquid height and increasing jet diameter for a given tank 
size. The average rms turbulent velocity is incorporated in 
Brown et al. (1990) condensation correlation to predict the 
condensation of vapor on a liquid surface. The results are 
in good agreement with available condensation data. 


1 . INTRODUCTION 

The interface condensation process which is principally con- 
trolled by liquid motion in s partially filled tank is of 
interest to space-based systems. In space environment, a 
preferred method to control cryogenic storage tank pressure 
is the use of jet-induced mixing (Aydelott, et al., 1985; 

Poth and Van Hook, 1972). For a laminar jet the interface 
condensation rate could be simply determined by the jet 
volume flow rate and jet subcooling (Lin, 1989; Hasan and 
Lin, 1990; and Lin and Hasan, 1990). However, the space- 
based systems are dimensionally large. The axial jet mixing 
system in most applications will be turbulent. Sonin et al. 
(1986); Brown et al (1989); and Brown et al . (1990) obtained 

an empirical correlation for the rate at which a pure vapor 
condenses on the free surface of a turbulent liquid in a 
steam-water system, for conditions where buoyancy effects are 
insignificant . 

Figure 1 shows the physical system and the* coordinates simi- 
lar to that used in Brown and Sonin's steam-water condensa- 
tion experiments (Sonin, et al., 1986; Brown, et al., 1989). 
The system consists of a vertical cylindrical tank of diame- 
ter D, partially filled with water. A steady turbulence is 
created in the water by an axial jer which is directed toward 
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FIGURE 1. Physical system and coordinates. 

the free surface from below, with the water circulating in a 
closed loop (not shown here) through a pump and cooler. 

Thus, a certain degree of subcooled axial jet flow can be 
maintained to induce vapor condensation at the surface. It 
was shown by Brown et al. (1990) that, for high jet submer- 
gence H/D > 3.1 and at conditions where buoyancy effect are 
insignificant, the condensation mass flux m of the pure 
vapor at the turbulent surface of its liquid 0 is given by 


St = 0.0198 


1 - — 

2 


Pr 


-0.33 


( 1 ) 


The condensation Stanton number is based on the condensation 
heat flux (m = h tg ), the local bulk liquid subcooling 
(T 3 “ T out)' ana the local rms value of turbulent veloc- 
ity (v s ' ) extrapolated from the bulk to the surface. For 
high jet submergence range 3.1 < H/D < 4.2, a correlation for 
the average value of rms turbulant velocity over the inter- 
face, v s ' , was given by Brown et al. ( 1990): 
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Brown et al. (1989) also conducted experiments to measure 
the mass condensation rate for jet submergence range 
0.5 < H/D < 3. However, the turbulence intensity was not 
measured. Following Thomas (1979), they modeled the interfa- 
cial turbulence intensity distribution and obtained an 
expression for the average value of v s ' over the interface 
as 


_ u d v H ( 3 ) 

TP = _1_ 10.4/3, - (7.14)5, - 3.06ft)- 1 ; 

where (3 and (3 are constants to be determined. Assuming 
Eq. (1) also applies to the low jet submergence, the conden- 
sation data of Brown et al. (1989) yield 

ft = 0.34 ft = 0.2 4 ( 4 ) 


The purpose of this study is to numerically determine the 

value of TP by exactly solving for the entire flow field 
in a tank with 0.5 < H/D < 3. The liquid-side turbulence as 
a function of relevant parameters is determined and incorpo- 
rated in Eq. (1) to predict the vapor condensation rate and 
compared with the experimental data of Brown et al. (1989). 


2 . PROBLEM FORMULATION AND MATHEMATICAL MODELING 

The physical system and coordinates are shown in Fig. 1. An 

axial liquid jet of diameter a and of velocity u., located 

at the tank bottom, is continuously introduced into the tank 
along the tank centerline. Liquid is withdrawn from the 
outer portion of the tank botrom with a withdrawing area of 
A . The volume flow rate of the liquid withdrawn is the 
same as that of the injected liquid jet such that the liquid 

level height H of the tank (or jet submergence) is main- 

tained constant. The effect of vapor motion on the liquid 
turbulence is assumed negligible. Also the liquid surface is 
assumed to be flat. 

The jet-induced mixing problem considered in the present 
study is steady-state and incompressible with gravity acting 
in the vertical negative-x direction. The dimensionless form 
of continuity and momentum equations with K-e turbulence 
model are expressed as 


3u* 1 3(rV) 

3x* r* 3r* 
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where the turbulence production term, G, is defined as 
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and the empirical constants associated with the turbulence 
model are 

C L = 1.44, C 2 = 1.92, and C u = 0.09 (11) 


Boundary conditions are required to solve the partial differ- 
ential Eqs . (5) to (9). At the centerline, the symmetric 

conditions are used. No-slip conditions for velocities are 
applied at the solid walls. Wall functions (Launder and 
Spalding, 1974) are used to calculate shear stress and are 
included as a source term in the governing equation for the 
velocity component parallel to the wall. For the axial jet, 
uniform velocity is assumed. The turbulent kinetic energy of 
the jet is set to be uniform at a specified level. The 
distribution of the dissipation rate is given by 

* 3/2 

* 3/4 K 3 (12) 

= C„ — 


where 1 is the mixing length scale (Escudier, 1966). 

CD 

Uniform velocity at the liquid-withdrawn plane is also 
assumed : 


u 


and v ouc = 0 


(13) 


At the free surface, zero axial velocity and shear free 
conditions are used: 


u = 0 and 

3 


3v * 
9x * 


= 0 


(14) 


Net fluxes of turbulent quantities K* and e* are assumed 
zero at the free surface and liquid withdrawn-plane: 
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The governing equations are numerically solved by a finite- 
difference method. Nonuniform grid distribution is used with 
concentration of the grid nodes in the centerline, near-wall, 
and near-interface regions where the gradients of flow prop- 
erties are expected to be large. The essential features of 
the numerical scheme are described in Lin (1989) and Lin and 
Hasan ( 1990 ) . 

The grid patterns of 60 by 41 and 84 by 41 have been shown to 
give reasonably grid-independent solutions for Ar < 1 and 
2.5 > Ar > 1, respectively. The convergent solutions are 
considered to be reached when the maximum of absolute resid- 
ual sums for u*, v*, and K* variables is less than 10' 6 . 

In the present study, calculations are performed to cover a 
wide range of the parameters Re., D/d and Ar: 10 000 
< Re^ < 60 000, 15 < D/d < 70, and 0.5 < H/D < 2.5. The 
value of K is fixed at 0.005 and A /A, = 80. 

J out' j 


3. RESULTS AND ANALYSIS 

The effects of the jet Reynolds number (Re ), aspect ratio 
(Ar) and tank to jet diameter ratio (D/d) on the turbulent 
quantities at and near the interface are of particular inter- 
est. Numerical solutions show that the dimensionless turbu- 
lent kinetic energy distribution (K*) at the interface is 
slightly increasing with Re j . However, this increase is so 

insignificant that the average value of K s * over the inter- 
face is essentially independent of Re. as 3 shown in Fig. 2. 
Figure 3 shows the distribution of dimensionless turbulent 

kinetic energy at the interface (K s ’) as a function of aspect 



Jet Reynolds number, Rej Radius, r" 


FIGURE 2. Average turbulent FIGURE 3. Effects of tank 

kinetic energy at the inter- aspect ratio and tank to jet 

face as a function of jet diameter ratio on the turbu- 

Reynolds number (H/D = 0.967, lent kinetic energy at the 
D/d = 23.6). interface (Re, = 20 000). 
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ratio, Ar, and tank to jet diameter ratio. D/d. An increase 
in either Ar or D/d results in a lower level of turbulent 
kinetic energy at the interface. This is because, for 
greater Ar or D/d, more energy is dissipated and less 
energy remains when the interface is reached. Also, greater 
values of Ar and D/d yield more uniform radial distribu- 
tion of K $ . For Ar greater than 3, the rms turbulent veloc- 
ities were found to be nearly uniform at and near the 
interface (Sonin et al., 1986). 

Komori et al. (1982) measured the turbulence intensities, 
i.e., rms values of three components of velocity fluctua- 
tions, in an open channel flow. It was observed that near 
the free surface, the intensity of vertical fluctuation 
decreases while the fluctuations parallel to the surface are 
promoted. This behavior has also been observed by Sonin 
et al . ( 1986 ) in a flow system similar to Fig. 1. At the 

interface, since the axial component of fluctuating velocity 
is damped out, only radial and circumferential turbulent 
velocities contribute to the turbulent kinetic energy. From 
the measurements of Komori et al. (1982), it can be seen that 
the turbulent intensities parallel to the free surface are 
nearly equal at the free surface. Thus, for an open channel 
flow, the rms value of the turbulent velocity at the free 
surface is approximately equal to the square-root of its 

turbulent kinetic energy. The average value of v' * over 
the interface can be calculated by 
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r 



*L/2 

K. r 
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(16) 


However, for the axial jet problem in the present study, the 
two components of interfacial turbulent velocity in the 
direction parallel to the interface may be significantly 
different. This observation is based on the experiments of a 
turbulent radial wall jet performed by Poreh et al. (1967). 
Therefore, the averaging process of Eq. (16) may not be 
appropriate. An alternative is to calculate the average 
value of the dimensionless turbulent kinetic energy over the 
interface by 

? = 8 f°- 5 K s V dr* (17) 

9 Jo 3 


and obtain the average value of rms dimensionless turbulent 
velocity by 


Equation (18) can be considered as a system value of rms 
dimensionless turbulent velocity at the interface and is used 


as the definition of v 


/ * 


in the present study. 
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It is 



noted that if the distribution of K s * is uniform, Eqs . (16) 
and (18) should give the same value of v' * . 


Figures 4 and 5 show that v g * is decreasing with increas- 
ing D/d and Ar, respectively. Based on the numerical 
solutions, a correlation is obtained for the average value of 

rms interfacial turbulent velocity (v'*): 


v' 


* 


3 


= (2.39 


- 0 . 645 Ar 



(19) 


for 0.5 < Ar £ 2.5 and 15 < D/d < 70. Within this range, 
v ' * is essentially a linear function of d/D and Ar. For 
a given tank size (D), the value of v g ' * is linearly in- 
creasing with increasing jet diameter (d) and decreasing 
liquid height (H). 


The value of v' ( = u. v' *) is the turbulent quantity 
required in Brown et al. (1990) condensation rate correla- 
tion. The average condensation mass flux at the surface can 
be obtained by integrating Eq. (1) and is expressed as: 


m h 


fg 


pC St (T - T ) v' 

I p C ' S OUt ' 3 


( 20 ) 




FIGURE 4. Effect of tank to 
jet diameter ratio on the aver- 
age rms turbulent velocity at 
the interface (Re. = 30 000, 

H/D = 0.9 67 ) . 


FIGURE 5. Effect of tank 
aspect ratio on the average 
rms turbulent velocity at 
the interface (Re. = 20 000, 
D/d = 2 3.6) . 3 
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Substituting Eqs . (1) and (19) into Eq. (20), the dimension- 
less system condensation rate (II) becomes 

II = 0 .0602 - 0.0163 - ( 21 ) 

D 

where n = HT h £g D d Pr^V^C^ (T s - T out )(l - Ja/2 ) ] as 

defined in Brown et al. (1989). Figure 6 shows good agree- 
ment between the measurement of Brown et al. (1989) and 
Eq. (21). Therefore, Eq. (21), based on the definition of 

Eq. (18) to predict v' * can be incorporated with Brown 

et al. (1990) condensation correlation for the condensation 
rate of a vapor with low jet submergence. Based on Eq. (19), 
the corresponding values of (3 and /? 2 for Eq. (3) are 
given by 

P x = 0.3 3 and p l = 0.2 3 ( 22 ) 


which is in very good agreement with Eq. (4) in the model of 
Brown et al. (1989). 


6x1 O' 2 



FIGURE 6. Comparison of Eq. (21) and 
Experimental data of Brown et al. (1989). 



4. CONCLUSIONS 


For an axial jet-induced mixing system as shown in 
Fig. 1, numerical calculations based on the assumptions 
and ranges of parameters considered lead to the following 
conclusions : 

o The radial distribution of turbulent kinetic energy 
at the interface (K g ) is more uniform for greater Ar and 
D/d with the maximum value located at the centerline. The 
average value of dimensionless turbulent kinetic energy over 
the interface is essentially independent of jet Reynolds 
number (Re ) and decreases with increasing aspect ratio (Ar) 
and tank to jet diameter ratio (D/d). However, experimental 
data are required for further quantitative evaluation. 

o The average (system) value of rms turbulent velocity 

(v' *) over the interface is linearly increasing with 
decreasing Ar and (D/d) for 0.5 _< Ar _< 2.5 and 15 < D/d 
< 70. Equation (19) for the prediction of v # ' can be 

incorporated with Brown et al. (1990) condensation correla- 
tion (Eq. (1)), for the vapor condensation rate in a tank 
mixed by an axial turbulent jet. 


NOMENCLATURE 

A W A surface area of the jet and outflow 

j ' out J 

Ar , B aspect ratio, H/D, tank to jet diameter ratio, 


C 

p 

d, D 

g 

H 

h, 

eg 

Ja 

* 

K 

1 

m 

m 

c 

Pr 


P'P g 

ft 

P 


Re 


i 


St 

c 


D/d 

specific heat at constant pressure 
jet diameter, tank diameter 
acceleration due to gravity 
jet submergence depth 
latent heat of condensation 
Jakob number, C (T -T , )/h. 

' p x 3 out 9 • ?g 

* / 2 

dimensionless turbulent K.E., R = K/u j 

mixing length scale 

local condensation mass flux 

Prandtl number of bulk liquid 

pressure, equilibrium hydrostatic pressure 

, , 2 

dimensionless pressure, p = (p - p )/(pu 3 ) 
jet volume flow rate 

dimensionless radial coordinate, r/D 
jet Reynolds number, pu d/p 

local condensation Stanton number, 
m h, /[pC ( T - T ) v ' ] 

c tq' 1 • p v s out' s J 
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V T out 


U 

V 

v 5 * 7 

v 7 ' 


,T temperature of the jet, outflow, and interface, 

respectively 

dimensionless axial mean velocity, u/u^ 
dimensionless radial mean velocity, v/u^ 
rms turbulent velocity 

dimensionless rms turbulent velocity, v > /u i 
dimensionless axial coordinate, x/D 


Greek symbols 


n 

P 

°k 

<7 

e 

Subscripts 

j 

out 

s 


dimensionless dissipation rate of turbulent 
kinetic energy, e/(uj/D) 

dimensionless system condensation rate 
dynamic viscosity 
empirical constants in Eq. (3) 
liquid density 

Prandtl number in K equation, cr R = 1.0 
Prandtl number in e equation, 


<7 =1.23 

e 


evaluated at jet inlet 
evaluated at outflow location 
evaluated at liquid-vapor interface 
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